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Xyt  general  analysis  is  presented  of  free  electron  lasers  in  which  a  static  periodic  magnetic  pump 
field  is  scattered  from  a  relativistic  electron  beam.  The  steady  state  formulation  of  the  problem 
is  fully  relativistic  and  contains  beam  thermal  effects.  Growth  rates  associated  with  the  radiation 
field,  efficiencies,  and  saturated  field  amplitudes  are  derived  for  various  modes  of  operation. 
Effects  of  space  charge  on  the  scattering  process  are  included  and  shown  to  play  a  dominant  role 
in  certain  situations.  Scaling  laws  for  the  growth  rates  and  efficiencies  at  a  fixed  radiation 
frequency  as  a  function  of  the  magnetic  pump  amplitude  are  obtained.  The  shear  in  beam  — ^  f- 
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The  Theory  of  Free  Electron  Lasers 


r.-  l  T 


I.  INTRODUCTION 


The  class  of  free  electron  lasers  (FELs)  in  which  a  pump  field  is  scattered  from  a  relativis- 

•  •  :  -i  ■'  •  .  '  ■  '  r  *  • 

tic  electron  beam  is  of  great  interest  as  a  potential  high-power,  tunable  source  of  coherent  radi¬ 
ation,  particularly  in  the  infrared,  visible  and  ultra-violet  spectral  regions.  The  concept  involves 

,!  ,  ..  «...  ,  4-  r,  ...  ;  »*  w  '  .  -  .  .  .  .  •.  I,  .  . 

the  stimulated  backscatter  of  a  pump  wave  from  a  relativistic  electron  beam.  The  pump  wave 
may  be  either  an  electromagnetic  wave  or  a  static  periodic  electric  or  magnetic  field.  For  a  stat¬ 
ic  periodic  pump  wave,  the  backscattered  radiation  frequency  from  a  relativistic  electron  beam 
is  a»  ~  (l+/3j)y^rc  *,  =s.2y,2c*»,  where  cp,  -  vr  is  the  axial  drifting  beam  velocity, 
Jz  ”  0 -p-)~'n,  k„  ~2w/l  and  ./  is.  the  period  of  the  pump, wave.  The  modulated  source 
currents  for  the  coherent  scattered  radiation  are  generated  ^  by  axial  bunching  of  the  electron 
beam  at  the  radiation  wavelength  through  coupling  of  the  scattered  waves  and  the  pump  field. 
The  mechanism  responsible  foj  this, axial  bunching,  is  the  ponderomotive  force  acting  on  the 
electrons  in  the  combined  fields  of  the  pump  and  radiation  waves. 


Analysis  and  design  considerations  pertaining  to  the  single  particle  scattering  process  have 
been  carried  out,  both  classically 1-14  and  quantum-mechanically. 15-17  When  the  electron  beam 
is  sufficiently  intense,  collective  effects  become  important  and  indeed  may  dominate  the  pro- 
cess.  Linear  analysis  of  the  FEL  have  been  performed  including  collective  effects,4*12, 14,1,-24 

*  .  *  '  .  *  •  '  i  ■’  {U  j  ,  »  .*  /  .r 

and  scattering  efficiencies  have  been  derived  for  various  FEL  scattering  regimes.5  5,12,  14,22 
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Free  electron  laser  experiments  with  pulsed  intense  relativistic  electron  beams  have  been 
conducted  at  a  number  of  laboratories.”-10  Submillimeter  radiation  at  MW  power  levels  were 
generated  with  electron  beams  of  energies  up  to  a  few  MeV  and  currents  in  the  multi  kA  range. 
In  these  experiments  collective  effects  play  a  dominant  role  in  the  scattering  mechanism  be¬ 
cause  of  the  high  beam  currents. 

In  another  class  of  FEL  experiments  at  Stanford  University,1112  relatively  low  current, 
high  energy  beams  were  employed  (/0— 1A,  24MeV  <  Eo  <  43MeV).  In  these  tenuous  beam 
experiments  collective  effects  are  negligible  and  single-particle  scattering  physics  apply.  Operat¬ 
ing  in  the  oscillator  mode,12  peak  powers  of  —7 kW  at  3.4  nm  were  generated  with  an  efficiency 
of  -  0.01%. 

In  this  paper  we  present  a  general  analysis  of  the  FEL  process  utilizing  a  right  handed  cir¬ 
cularly  polarized,  spatially  periodic  magnetic  pump.  A  schematic  of  the  FEL  configuration  is 
shown  in  Fig.  (1).  The  analysis  is  fully  relativistic  and  performed  explicitly  in  the  laboratory 
frame.  Our  formulation  shows  that  depending  on  the  beam  and  pump  parameters,  several  dis¬ 
tinct  interaction  processes  can  be  distinguished.  Our  results  are  applicable  to  both  the  tenuous 
and  intense  beam  type  experiments.  Growth  rates  (or  gains)  together  with  saturation 
efficiencies  are  derived  for  the  various  FEL  regimes.  A  condition  for  the  neglect  of  collective 
effects  for  the  low  gain  FEL  process  is  derived.  Scaling  laws  for  the  growth  rates  and 
efficiencies  at  a  fixed  output  frequency,  as  a  function  of  the  pump  amplitude  are  given.  The  de- 

«.  .  i 

trimental  effect  of  axial  velocity  shear  on  the  beam  due  to  self  fields  is  discussed  and  various 
methods  of  reducing  this  shear  are  suggested.  In  addition,  an  illustration  of  a  far  infrared  two- 
stage  FEL  using  a  3  MeV  electron  beam  is  presented.  Here  the  output  radiation  wavelength  is 
decreased  approximately  by  the  factor  tyf  compared  to  the  pump  wavelength  instead  of  the 
factor  2  y  ?  for  a  single  stage  FEL. 
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II.  DERIVATION  OF  GENERAL  FEL  DISPERSION  RELATION 


The  pump  is  chosen  to  be  a  right  handed  circularly  polarized  magnetic  field  given  for  z  > 

0  by 

B„  -  B„(excos{k„z)  +  e„sin0c„z)),  (1) 

where  B„  is  constant  and  k„  -  2 n/l,  see  Fig.  (1).  The  representation  of  the  pump  field  in  (1) 
is  a  good  approximation  near  the  axis  of  an  appropriate  coil  winding.  The  vector  potential  asso¬ 
ciated  with  B„  is  A„  -  -B Jk„.  For  particles  in  the  field  given  by  (1),  the  canonical  momenta 
in  the  x  and  y  directions  as  well  as  the  total  momentum  are  constants  of  the  motion  and  are 
given  respectively  by 


a(px,  z)  -  px  -  A,a(z).  (2a) 

c 

P(p„.  z)  -  Py  -  A,„(z),  (2b) 

c 


h(p)  -  Ip!. 


(2c) 


where  p  is  the  momentum. 

We  assume  that  the  interaction  between  the  relativistic  electron  beam  and  the  pump  field 
has  reached  the  temporal  steady  state  so  that  the  radiation  fields  are  proportional  to  exp(-rW) 
where  w  is  the  frequency  of  radiation.  The  radiation  and  space  charge  fields  are  given  by 


E(z,  ,)-A 


f:  + if- A  (z)  (ex  +  /pv)|p—'  +  c.c. 


B(z.  () 


/  9A  (z) 
2  9z 


(f,  A-  /e,)p-'-'  +  c.c. 


(3a,b) 
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where  the  associated  potentials  are  *(z./)  -  0/2)  *(z>  exp(-/«/)  +  c.c.  and  A(z.r)  -  (1/2) 
A  (z'  (ex  +  iey)  exp i-imt)  +  c.c.. 

Using  the  relativistic  Vlasov  equation  we  expand  the  electron  distribution  function  to  first 
order  in  the  scattered  fields  E  and  B  about  its  equilibrium.  It  proves  very  convenient  first  to 
transform  the  independent  variable  (p,  z,  t)  of  the  distribution  function  to  the  new  indepen¬ 
dent  variable  (o.  p.  u.  z,  t).  The  electron  distribution  is  then  written  as 

g(a.  p.  u.  z.  I )  -  gm(a.  P.  u )  +  p,  u,  z,  l)  (4) 

where  g(m(a.  p.  u)  is  the  equilibrium  distribution  and  is  an  arbitrary  function  of  the  constants 
of  the  motion  and  gu)(a,  p,  u,  z,  I )  is  the  perturbed  part  of  the  distribution  which  is  propor¬ 
tional  to  either  4>  or  A.  It  is  straightforward  to  show  that  the  perturbed  part  of  the  Vlasov 
equation  for  g(l>  is 


where 


dg(l)  +  P:  dgIU 
dt  ym„  dz 


-  H  e~"“  +  c.c. 


(5) 


H  -  H(a.  P.  u.  z) 


+ 


A 


u 


(6) 


and  the  dependent  variables  are 


v.  -  v.(a,  P,  u,  z)  -  P_-(a.  P.  u.  z)Hy(u)  m„). 


P- 


-  pAa,  P,  z)  -  Px  +  iPv  -  a  +  >P  ~ 


I *\B. 

ck„ 
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Pz  -  pz(a,  P.  U,  z)  -  ( U 2  -  p}  -  p}V'2, 

\e\ 

Px  "  Px(«.  z)  «  a  +  -L^L  ^ox(z), 

Py  "  Py(P.  Z)  “  P  +  ^oy(z), 

y  -  y(«)  -  (1  +  u2/(m2c2))  'n.  <7a-0 

The  general  solution  of  Eq.  (S)  is 


*(,)  .  !<■>  e-/-»  + 


where 


#<u  -  #<l)(o,  P.  u.  z)  -  f  m  dz'M(a,  p.  u,  z,  z')  H(a,  p,  u,  z'), 


(8) 


M(a,  P,  u.  z,  z')  - 


giiria.fi.u.t.z') 

Vj(a,  P,  u.  z') 


and 


r(«,  p,  u,  z,  z') 


dz" 

vt(a,  P,  u.  z')  ’ 


The  beam  has  been  taken  to  be  unperturbed  at  z  -  i.e.,  g<n  (a,  p.  u,  -°°)  -  0. 

Rearranging  Eq.  (8),  the  Fourier  transform  of  the  perturbed  part  of  the  distribution  func¬ 
tion  becomes 


#<u(a,  p.  u,  z)  -  |</+(ar, 0,w,z)  +  /^-J  +  6,(0.  P,  u,  z)~J#<0)(a,  fi,  u ),  (9) 


where 


*m>vt  ,napRisMn|r 
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6+  -  M  J_m  dz'  M(a,  P,  u,  2,  z')  J/  w  -  Vj(a,  P,  u,  z‘)  ~J  A(z‘) 


—Mi 


A(z). 


6,  -  M  —■  f  m  dz’  M(a,  P,  u,  z.  z')  {/«/?_(«,  p,  z')  A  (z') 


-cpz(a,  P,  u,  z')  -  . 


In  Eq.  (10a)  we  have  integrated  by  parts,  using  the  fact  that  the  radiation  field  vanishes  at 
z  “  — °°.  The  expression  for  #n>  in  Eq.  (9)  determines  the  first  order  perturbation  of  the  elec¬ 
tron  distribution  function  due  to  scattered  fields  with  arbitrary  axial  spatial  dependence,  is 
correct  to  all  orders  in  the  pump  field  amplitude,  and  contains  thermal  effects  in  g„. 


The  perturbed  current  density  which  drives  the  scattered  fields  is  given  by 


J{z,  t)  —  (J+(z)  e+  +  Jz(z)  +  c.c. 


where 


f- f-  ['♦<«•  >■  ’>  I  (12) 

):{z)  m„  Jo  j—  j—  p  |p..(a,  p,  u,  r)J  y(u)p;(a,  p,  u,  z)‘ 

•  ;  •  *  '  J  ’  » 

where  />+(«.  p,  z)  -  px  -  ipy  -  a  -  ip  -  (|*|  Bjck„)  exp(-/*„r),  and  e+  s  (c,  +  *,)/!. 


To  evaluate  the  current  density  in  (12),  we  take 


g(0,(«.  P.  u)  -  n„6(a)  t(p)  g.(u). 


(13) 
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where  n„  is  the  unperturbed  beam  density,  assumed  to  be  uniform  in  space,  and  ga(u)  is  arbi¬ 
trary  but  subject  to  the  normalization  condition  J"o  du  g„(u)  u/uz  *■  1.  The  delta  functions  for 
a  and  /9  arise  from  assuming  that  the  equilibrium  transverse  momentum  is  due  solely  to  the 
pump  field,  i.e.,  that  transverse  thermal  effects  can  be  neglected.  Using  the  distribution  func¬ 
tion  (13),  we  find,  after  some  lengthy  algebra,  that  the  Fourier  coefficients  of  the  current  densi¬ 
ties  given  in  (12)  are 


<»h  r" 

4  it  Jo 


—  (1  -  ^-)  A(z) 
yuzc  2 


\ 

+  +  A  (14., 


7  l.\  _  °>h  —  r~  J..  J f 2  J_i  d$(z')  —ivz'/v, 

J-iz)  8^  m"  Jo  due  Jo  dz  —ftT  e  +  ~T~ 


—  A  (z')  e~'*M  dgjdu. 
C 


where  -  (4ir«„  |e|J/m„)l/2,  /3 x  -  ps(u)  -  ft„/(y(u)  v:(u)k„),  Sl„  -  \e\BJm„c, 

*P(z)  -  («/ v.(u)  -  k„)z,  u.  -  u.(u)  -  p.(a  -  0,  0  -  0,  u.  z)  -  ( u 2  -  m„2n  2//fc02)l''J.  and 

“  u:/(y(u)m„).  The  limits  of  integration  over  z'are  from  0  to  z  and  not  from  —  oo  to  z,  be¬ 
cause  the  amplitude  of  the  various  fields,  i.e.,  A„,  <j>  and  A,  are  assumed  to  build  up  from  zero 
at  z  —  -oo  to  their  initial  values  at  z  —  0  in  a  distance  which  is  small  compared  to 
(*  +  +  k„  -  «/v-)-\  where  *+  is  the  wavenumber  associated  with  A(z).  The  limits  of  integra¬ 
tion  over  z'can  therefore  be  changed  from  (— oo,  *)  to  (0,  z)  without  loss  of  accuracy.  Because 
the  characteristic  length  (k  +  k„  -  at/v.)~l  is  much  longer  than  the  wavelength  of  the  pump 
field,  this  situation  is  necessarily  satisfied  in  any  experimental  configuration.  The  driving 
current  density  in  Eqs  (14)  contain:  (i)  the  ponderomotive  potential  manifested  in  flie  term 
PiA(.z)\  (ii)  collective  effects  from  the  scalar  potential;  (iii)  arbitrary  axial  variation  of  the 
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excited  fields  0  and  A;  (iv)  ballistic  terms  propagating  from  the  boundary  at  z  —  0  associated 
with  the  lower  limit  on  the  z'  integral;  and  (v)  arbitrary  thermal  nature  or  the  beam  manifested 
in  g„(u). 

The  analysis  is  closed  by  taking  the  perturbed  current  density  of  Eqs.  (14)  to  be  the 
source  current  in  the  wave  equations  for  A  and  <f>: 

1-^-  +  'A  (z>  “  —  J+<*>;  (15a) 

[or  c  c 

^  osb) 

or  ai 

A  number  of  different  scattering  regimes  can  be  distinguished  using  the  general  form  for  the 
driving  currents  expressed  in  Eqs.  (14).  We  shall  discuss  in  detail  those  regimes  which  appear 
to  be  important  for  the  development  of  efficient,  high-power  FEL's. 

III.  LOW  GAIN  LIMIT 

The  first  case  we  consider  is  the  low  gain  or  short  cavity  regime,  where  collective  effects 
do  not  play  a  dominant  role  and  the  electromagnetic  wave  is  nearly  of  constant  amplitude.  By 
low  gain  limit  we  mean  that  the  total  integrated  gain  of  the  radiation  field  is  much  less  than  un¬ 
ity.  This  limit  corresponds  to  the  parameter  regime  of  the  experiments  carried  out  at  Stanford 
University51  52  with  highly  relativistic  (<48  MeV),  low  current  (<  2.4 A)  beams.  Neglecting 
collective  effects  implies  that  0(z)  «  A  (z)  </3I>;  the  condition  on  the  beam  density  for  this 
inequality  to  be  satisfied  is  given  below.  Taking  the  electromagnetic  field  to  be  of  the  form 

A(z)-A(  O)e'fok*(:),k  (16) 

where  |lm(fc+)|  <<  R e(*+)  and  expj^  <fe'Im(fc+(z'))  l.  With  this  representation,  together 
with  (14a)  and  (15a),  the  dispersion  relation  takes  the  form 
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^o+  -  «*2/c2  +  2<c0+  8A(z) 


-oil 


r* 

Jo 


u  ,  mn  n-*,w^(*o+  +  *o>K  ft 

—  (l+tf/2)  -  -^W  -(1  -g  „ - — 1  A 

•y«2  2  ot-vz(k0+  +  k0)  du 


g0(u),  (17) 


where  /c+(z)  —  A0+  +  8/c(z)  ,  |8<c|  «  k0+  and  k0+  =2 y?ck0  real  and  constant.  Solving  for 
the  imaginary  part  of  8A(z)  we  obtain 


ul/c2 
2k 


Im(8/c(z))  =  — 


o+ 


sin(<u/  —  (k0++kp))z  dg0 
w/vz-(k0+k0)  du  ' 


The  total  integrated  gain  in  the  wave  amplitude  over  the  interaction  region  of  length  L  is 
defined  as  GL  -  -  fQ  lm(8k(z'))e/z’.  It  is  straightforward  to  show  that  if  the  thermal  energy 
spread  of  the  beam,  is  such  that 

EJE0  <  <  y?0  ML  ,  (19) 


the  beam  can  be  considered  to  be  mono-energetic  in  Eq.  (18).  In  the  above  inequality 
£o  ”  (yo- 1)  mo*-2  =  ToW0c2  is  the  total  kinetic  energy  of  each  beam  particle, 
y.-o  “  (l-v-o/c2)-l/2  and  A  -  2n/k0+  is  the  wavelength  of  the  radiation  field.  For  future  use 
we  note  that  -  y0y;2o  mo  v:oKii  where  F„,  is  the  beam  thermal  velocity.  Assuming  (19)  to 
be  satisfied  we  can  use  the  distribution  function  g0(u)  -  ( u./u)8(u-u0 )  and  find  that  GL  is 


G,  - 


Q>ft/c2 

8<c0+ 


m„oiL2  dufil 


sin(<c0+-Kc0-(i>/v.-)£/2  2  dg0 


(.k0++k0—o)/vz)L/2  J  du 


sinf»n 


«  1 


(20) 


where  (  -  o»h/(y/y^c  k0),  y0  -  y(w0h  y:o  -  y:(«o).  /3ql  ”  &(«<>).  -  (w/v; o  -  *0+-*„)Z./2, 

and  u o  is  the  magnitude  of  the  total  particle  momentum.  The  function  8(sin0o /0„)J/80O  has  a 
maximum  value  of  0.54  when  0O-  —1.3,  hence,  the  maximum  total  gain  is 
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GLm„~((/*)Wk0L)3  (21) 

and  is  much  less  than  unity.  Using  Eq.  (ISb)  in  conjunction  with  Eq.  (14b)  we  And  that  the 
condition  for  neglecting  collective  effects,  i.e.,  4>(z)  «  A  (z)  <0±>  ,  can  be  written  as 


2w 


|  rm  .  (1-e  » _ )  too 

m°Uo  u  (o/vl  -  (*0+  +  <to)  3“ 


«  1. 


which  reduces  to 


*  L  k0  1 

2>io 


«  1. 


(22) 


The  L 2  dependence  in  (22)  is  due  to  the  dependence  of  the  perturbed  density  on  length.  The 
density  modulations  on  the  beam  can  be  shown  to  increase  as  z2  in  our  present  limit. 


It  will  be  necessary  to  obtain  the  difference  between  the  phase  velocity  of  the  ponderomo* 
live  wave  and  the  axial  electron  velocity  when  deriving  the  saturation  efficiency  in  Section  V. 
From  the  definition  of  Oq,  this  velocity  difference  is  simply 


V-  “  *'.-o  -  - 


<o 


*0+  +  *0 


90c 


yfoL  *o 


where  k0+  =  2y}0  k0  »  k0  and  the  value  of  0O  extends  from  0  to  = 
maximum  positive  gain. 


(23) 

-3  for  the  domain  of 


IV.  HIGH  GAIN  LIMIT 

In  contrast  to  the  first  case  we  now  consider  the  long  cavity  limit  where  the  excited  Aeld 
amplitude  spatially  exponentiates  several  times  within  the  interaction  region.  Under  these  con¬ 
ditions  the  terms  containing  the  boundary  conditions  at  z  «  0  can  be  neglected,  i.e.. 
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^ (0)  «  <£(2)  and  A  (0)  «  A  (z).  We  assume  that  conditions  appropriate  to  a  cold  beam  are 
satisfied,  hence  gu(u)  **  (u2/u)  8(«  -  «„);  we  shall  return  to  this  point  later.  The  potentials 
can  be  represented  by 


0(r)  -  *(0) 

(24a) 

A(z)  -A(0)e'*+z 

(24b) 

where  k  —  k+  +  ka  and  k+  is  complex.  Substituting  the  potentials  represented  by  (24)  into  Eqs 

(14) ,  in  conjunction  with  the  cold  beam  assumption,  and  making  use  of  the  wave  equations 

(15)  we  obtain  the  following  dispersion  relation  for  <■>  and  fc. 


1  2 

D(ot,  k  -  k0)  ((«  -  v„k)2  -  <o}/(y20y0))  -  — /)(«,  *) 

2y*  ck„ 


(25) 


where  D(w.  k)  -  co2  -  c2k2  -  a>i/y0 ,  y0  -  y(u0),  vzo  -  vz(u0)  and  yto  -  y2(u„).  The  elec¬ 
tromagnetic  wave  approximately  satisfies  the  dispersion  relation  D(io,  k+  -  k  -  k0)  =  0, 
hence,  we  can  replace  D{ w,  k)  on  the  right  hand  side  of  (25)  by  -2 kk0c2.  Also,  since 
*+  =  <o/c  we  approximate  D(o»,  k+)  by  -2 c2k+{k+  -  («J  -  <o}/y B)lll/c).  The  dispersion  rela¬ 
tion  can  now  be  put  into  the  simple  form 


(k  -  (K  +  k„ ))  ( k  -  (w/v„,  +  k))  ( k  -  («/vro  -  k ))  ~  — y  k0 


(26) 


where 


K  -  (<u2  -  o>l/y0)U2/c, 

*  “  tob/bzo  yto  yy2), 
a2  -  ((l„/ck„)2  (wj/(y2  v£))  -  (f0ai*o)2, 
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001  "  ^  o/^y  0^20^0^  • 

and 

(  “  cka) . 

Further  simplification  of  (26)  is  obtained  by  setting 

k  -  tt/vzo  +  a  +  8 k  (27) 

where  Bk  is  in  general  complex  and  |8  k\  «  k.  Since  vm  *=  cand  <u  »  «*/>/r <J,  we  find  that 
K  =  oi/c  -  (2k0/(4yl.)t  k  =  £*„/>„  so  that  (26)  reduces  to 

8* (8*  +  llkjyj  (6k  -  A *)  -  -a1*,#  (28) 

where  A *-•*«-  u/(2cyzo)  is  a  mismatch  parameter.  At  this  point  it  is  convenient  to  evaluate 
the  difference  between  the  phase  velocity  of  the  longitudinal  wave  and  the  initial  axial  beam 
velocity  in  the  high  gain  limit.  This  velocity  difference  is  given  by 

V  ~  vtn  =  -Av  -  WReOc)  -  v.„ 

_  -(k  -f  Re(8lc))  c  ^9) 

:n  k„ 

where  we  have  used  the  expression  for  k  in  (27).  The  expression  for  Av  in  (29)  will  be  used 
later  to  obtain  an  estimate  for  the  saturation  efficiency  and  maximum  radiation  field.  We  now 
discern  two  important  limits  of  the  dispersion  relation  (8). 

a)  Weak  Pump  Limit 

For  a  pump  magnetic  fteld  amplitude  such  that  «  £r„,  s  4((/y*,)'n  the  space 
charge  potential  dominates  the  ponderomotive  potential  and  collective  effects  are  important.14 
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That  is,  in  this  case  the  electron  susceptibility  \  “  —  (<o  —  v.-„k)~2w2/(y2,y,;)  is  approximately 
— 1  and  the  electrostatic  wave  is  nearly  an  eigenmode  of  the  system.  This  regime  of  scattering 
corresponds  to  setting  |SAc  |  «  2 £kjy.„  in  the  dispersion  relation  (28),  which  becomes 

6k  (8k  -  A*)  «  -a2y:„/(40  (30) 


with  the  growing  root  given  by 


(31) 


The  condition  for  instability  is  clearly  a2y.„/£  >  (A k)2  and  the  maximum  spatial  growth  rate 
occurs  when  A/r  -  0  and  is 


fin»«  —  “  2  y/£y ;ti  k„  (32) 

Using  (32)  we  see  that  the  condition  |8/c|  «  2 (k„/y:o  is  equivalent  to  the  weak  pump  condi¬ 
tion,  i.e.,  /Soa  «  fieri,.  In  this  FEL  regime  we  find  that 


(tkjyz„  +  AA/2) 

vph  -  vzo  -  -Av-  - - -  — - c 

2y:„k„ 


where  (31)  has  been  used  for  Re(8*),  and  A/c  <  0<u*«- 


b)  Strong  Pump  Limit 


(33) 


In  this  regime  the  pump  magnetic  field  amplitude  is  sufficiently  strong  to  satisfy  the  in¬ 
equality  fi„x  »  ficri,  =4((/y i,)U2.  The  ponderomotive  potential,  which  is  proportional  to  the 
pump  amplitude,  completely  dominates  the  space  charge  potential  in  the  strong  pump  regime14 
and  |xl  «  1.  This  is  a  single  particle  scattering  regime  where  collective  effects  are  negligible. 
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For  p„x  »  Pml  we  neglect  2(kjy:„  compared  with  8A  in  (28)  and  the  dispersion  relation  be¬ 
comes 

(8 A)2  (8A  -  AA)  -  —a2k„/2.  (34) 

The  maximum  spatial  linear  growth  rate  occurs  Tor  exact  Frequency  matching,  i.e.,  A k  —  0,  and 
is  given  by 

I™,  s  -Imtt*)™.,  -  ij^M^'k,,.  (35) 

while  at  this  Frequency  Re(8A)  —  (£/9,a)2/3A„/24/3.  The  real  part  oF  6k,  which  is  a  Function  oF 
A  k,  has  the  maximum  value  Re(8A)mil,  -  (f/3,a)J/3A„  when  A  *  -  y«0„i)2/3A„.  The  velocity 

difference  in  (29)  also  attains  its  maximum  value  when  Re (8 A)  -  Re(8A)m„  which  is  the  point 
where  the  growth  rate  vanishes.  As  we  shall  see  in  the  next  section,  the  energy  extraction  is 
proportional  to  Av,  and  hence  the  maximum  efficiency  is  attained  close  to  the  point  oF  vanish¬ 
ing  growth  rate. 

V.  SATURATION  LEVELS 

To  obtain  estimates  For  the  saturation  levels  in  the  different  FEL  regimes  we  resort  to  ar¬ 
guments  based  on  electron  trapping.14  In  the  cold  beam  limit,  we  assume  that  saturation  occurs 
when  the  beam  electrons  become  trapped  in  the  total  longitudinal  wave,  i.e.,  space  charge  plus 
ponderomotive  potential.  The  difference  between  the  longitudinal  wave  phase  velocity  and  the 
axial  electron  velocity  is  initially  v^,  -  v„,  -  -Av,  where  the  difference  Av  is  greater  than  zero 
For  instability  and  depends  on  the  particular  FEL  regime  as  well  the  Frequency  mismatch  param¬ 
eter  AA  (see  Eq.  (32)].  Assuming  all  the  particles  to  be  deeply  trapped,  we 


14 
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may  estimate  that  at  saturation  vp/l  -  v!MI  ~  Av  where  vzsa,  is  the  average  axial  electron  veloci¬ 
ty  at  saturation  and  vp/l  is  assumed  to  remain  fixed.  The  maximum  decrease  in  the  axial  beam 
velocity  is  ~2  A  v,  corresponding  to  a  change  of  particle  kinetic  energy  by  an  amount 

*Eke  =  (2A  v)m„c2 

”  -2y„  yl,  m„  v.n  A  v.  (36) 

The  energy  conversion  efficiency  is,  therefore, 

V  “  - - —  ~  ~  2y  }„  A  v/c.  (37) 

(y„  -  l)m0c3 

Similar  arguments  have  been  used  to  obtain  good  estimates  for  efficiency  in  two-steam  interac¬ 
tion  processes.33  The  vector  potential  at  saturation,  z  -  z„„  can  be  found  by  applying  the  con¬ 
servation  law  for  total  energy  flux.  The  result  is 

\A(z  -  zM,)|  -  |M*(z  -  0)|>  +  ^  -^]2  .  (38) 

In  the  low  gain  FEL  limit,  described  in  section  III,  the  efficiency  which  is  given  by  Eq. 
(37)  together  with  Eq.  (23)  is 


,  ' TT„ •  (39) 

for  the  highest  gain  band,  9„  ranges  from  0  to  -3.  The  maximum  gain  occurs  when  0„  —  -1.3 
and  is  given  in  Eq.  (21).  The  amplitude  of  the  vector  potential  at  saturation  is 

\A(z  -  z„,  -  £.)(  =  \A(z  -  0)|  (1  +  CE).  (40) 

Comparing  (38)  with  (40)  we  find  that  the  input  signal  needed  to  cause  saturation  at  z  «  L  is 


M,o)l“fe  t£|  ViF- 


with  Gt  given  by  Eq.  (21). 
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In  the  high  gain,  weak  pump  FEL  regime  discussed  in  section  IVa,  the  efficiency  is 

V  -  ((/Vza  +  A*/(2*„))  (42) 

where  we  have  used  Eq.  (33)  in  conjunction  with  (37),  noting  that  the.  mismatch  parameter  is 

A*  <  P  Equation  (42)  is  valid  in  the  high  gain,  weak  pump  parameter  regime  and 
hence,  the  second  term  is  somewhat  smaller  than  the  first.  The  amplitude  of  the  vector  poten¬ 
tial  at  saturation  in  this  case  is 

]A  (z " Zt0,)  1  ®  ( 2^t  TT )  v'/2,  (43) 

where  Eq.  (38)  was  used  together  with  the  condition  |>4  (z  -  zM,)|  »  |  (0)  | . 

Finally,  we  consider  the  high  gain-strong  pump  case.  The  efficiency,  using  Eqs.  (29)  and 
(37),  is  given  by 


v  "  p”1  +  ^1*  (44) 

where  Re(8k)  <  ((Pnj.)2/3k„,  the  equality  holding  where  the  spatial  growth  rate  vanishes. 
When  the  growth  rate  is  maximum  (see  Eq.  (33))  the  efficiency  is 


V  -  2-4/3({0„I)1/3  +  (/y„  (45) 

The  second  term  in  (44)  is  small  compared  to  the  first  in  the  strong  pump  limit.  The  saturated 

value  of  the  vector  potential  is  given  by  Eq.  (43)  together  with  (43). 

•  i 

VI.  GROWTH  RATE  (GAIN)  VS.  B0  FOR  FIXED  OUTPUT  FREQUENCY 

It  is  of  interest  to  determine  the  scaling  laws  for  the  growth  rate  (or  total  gain)  and 
efficiency,  for  a  fixed  output  frequency,  as  a  function  of  the  magnetic  pump  amplitude  B„.  To 
obtain  these  scaling  laws  for  a  Axed  output  frequency  «  -  2 ylck„,  i.e.,  fixed  y„  and  we 
note  that  the  total  gamma  can  be  written  as 
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y„  **  r-,,0  +  (|r|ff„/»i„r2A„)J),/i  (46) 

Therefore,  when  /)„  approaches  and  becomes  larger  than  the  critical  magnetic  field  amplitude, 

B,nl  -  (w„r2/|c|)  A„  -=  (10. 6// (on))  AC  (47) 

the  axial  gamma,  y.„  becomes  significantly  smaller  than  the  total  gamma  y„.  For  fixed  y.„  and 

k„  when  B„  «  Bc,„  the  total  gamma  is  nearly  equal  to  y.„;  when  B„  »  however, 

y  n  y  :n  B,,/ Biril . 

In  the  low  gain  FEL  regime,  Eq.  (21)  shows  that  the  maximum  total  gain  is  proportional 
to  (or  B„  «  B,n,  and  falls  off  as  B„  1  lor  B„  »  B,.r„.  The  efficiency  given  in  Eq.  (39)  is 
independent  of  the  pump  magnetic  field  amplitude  for  fixed  y  and  k„. 

In  the  high  gain-weak  pump  FEL  case  the  maximum  spatial  growth  rate,  see  Eq.  (32),  is 
proportional  to  B„  for  B„  «  Bir„  and  decreases  as  B~u*  for  B„  »  B,,,,.  The  efficiency  on 
the  other  hand,  see  Eq.  (42)  is  independent  or  B„  for  B„  «  B,r„  and  falls  off  as  B„  for 
B„  »  B,t„.  For  the  high  gain-strong  pump  case  Eqs.  (35)  and  (45)  shows  that  for  B„  «  Binl 
both  the  maximum  growth  rate  and  efficiency  increase  as  B,}n  whereas  for  B„  »  Btr,,  both  the 
growth  rate  and  efficiency  fall  off  as  B~'n. 

These  scaling  laws  for  fixed  output  frequency  indicate  that  for  all  the  FEL  regimes  which 
have  been  considered,  the  optimal  magnetic  pump  amplitude  is  one  where  B„  =  B,nl. 

VII.  DISCUSSION 

(a)  Energy  Shear 

In  the  preceding  formulation  of  FELs,  we  have  neglected  any  effects  or  energy  shear 
across  the  beam.  Such  a  shear  arises  owing  to  the  self  electrostatic  potential  drop  within  the 
beam.  This  leads  to  a  radial  dependence  of  the  beam  kinetic  energy  in  the  equilibrium  state. 
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The  energy  shear  results  in  a  shear  in  the  axial  equilibrium  velocity  across  the  beam  and,  there¬ 
fore,  is  equivalent  to  a  beam  temperature.  For  an  axially  propagating  beam  of  radius  r„  the 
effective  beam  temperature  is  of  order 

A£=s|e|A*- j-^j— j  £„ 

where  A4  is  the  self  potential  drop  across  the  beam  from  r  —  0  to  r  -  r„  and  £„ 

is  the  kinetic  energy  of  the  electrons.  A  necessary  condition  for  the  validity  of  the  cold  beam 

approximation,  in  all  the  FEL  regimes  which  have  been  considered,  is 

A £/£„  «  t>. 

This  inequality  may  be  invalid  at  sufficiently  high  beam  densities. 

A  more  refined  analysis  taking  account  of  the  energy  shear  should  also  consider  the  radial 
gradient  of  the  pump  field,  which  is  necessary  to  satisfy  0.  The  radial 

dependence  of  the  pump  produces  a  shear  in  the  eqiulibrium  transverse  velocity  which  will  tend 
to  compensate  for  the  shear  in  axial  velocity  due  to  self  field  effects.  Other  possible  approaches 
which  may  be  considered  to  eliminate  the  axial  velocity  shear  include  (i)  establishing  Brillouin 
flow  in  the  beam  by  applying  an  axial  magnetic  field  or  (ii)  creating  the  beam  on  a  non- 
equipotential  surface  so  that  the  applied  potential  shear  just  cancels  out  the  self  potential  shear. 
In  the  following  example  self-field  effects  will  be  neglected. 

(b)  Two  Stage  FEL 

As  an  illustration  of  a  far  infrared  radiation  source  we  consider  a  two-stage  FEL  generator. 
In  a  two-stage  FEL,  two  consecutive  and  distinct  scattering  processes  lake  place  within  a  single 
electron  beam.  The  output  radiation  from  the  first  stage,  in  which  the  pump  is  a  circularly  po¬ 
larized  static  magnetic  field,  is  reflected  back  on  the  beam  and  used  as  the  pump  wave  in  the 
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second  stage.  This  configuration  is  schematically  depicted  in  Fig.  (2).  The  -final  wavelength  of 
the  output  radiation,  from  the  second  stage  is  A  =  //8-y4  instead  of  Illy}  as  would  be  the  case 
in  a  single  stage  device.  Hence,  in  a  two  stage  FEL,  far  shorter  output  wavelengths  can  be  real¬ 
ized  for  the  same  electron  beam  energy.  The  pump  field  in  the  second  stage  is  a  circularly  po¬ 
larized  electromagnetic  wave  and  not  a  circularly  polarized  static  magnetic  field  as  in  Eq.  (I). 
Our  results  for  a  magnetic  pump,  however,  also  apply  to  a  circularly  polarizated  electromagnetic 
pump  if  the  electron  beam  is  highly  relativistic.  To  see  this  we  note  that  in  the  beam  frame, 
the  two  pump  waves  are  equivalent  if  we  set  B()2  -  2£,(r  -  z„„)  and  k02  -  2  X+,  where  B0 2  and 
k02  is  the  magnetic  field  amplitude  and  wavenumber  of  the  equivalent  magnetic  pump  in  the 
second  stage  and  E\(z  —  zw„)  and  X+)  is  the  saturated  electric  field  amplitude  and  wavenumber 
of  the  reflected  output  radiation  from  the  first  stage.  The  relevant  parameters  for  this  example 
are  contained  in  Table  I.  The  results  outlined  in  Table  I  demonstrate  that  in  principle  a  rather 
low  energy  electron  beam  (£„  —  3MeV)  is  necessary  to  generate  far  infrared  radiation  using  a  2 
cm  wavelength  magnetic  pump.  The  radiation  to  beam  power  efficiency  of  0.085%  may  be 
greatly  enhanced  by  adiabatically  varying  the  longitudinal  wavelength  of  the  electromagnetic 
pump  in  the  second  stage.  Contouring  the  pump  period  for  the  purpose  of  enhancing  efficiency 
has  been  suggested  in  Ref.  (14). 

Recent  non-linear  calculations  have  shown  that  efficiency  enhancement  factors  greater 
than  100  can  be  achieved  by  varying  the  wavelength  of  the  static  magnetic  pump  field.34  In  the 
case  of  an  electromagnetic  pump  the  axial  wavelength  may  be  contoured  by  varying  the 
waveguide  wall  radius.  Work  is  now  underway  at  the  Naval  Research  Laboratory  to  fully  evalu¬ 
ate  this  approach. 
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Table  I  —  Illustration  of  a  Far-Infrared  (X  -  2m  m) 
Two  Stage  FEL 


Electron  Beam 


Energy:  £0  «*  3  MeV  (y0  —  7),  Current:  /„  —  10  kA,  Radius:  r0  “  0.3  cm. 


First  Stage 

Second  Stage 

Pump  Amplitude 

Bor  -  5  kG 

B0 2  “  2£,(r  -  2ml)  -  15.7  kG 

Pump  Wavelength 

l0 1  —  2  cm 

I02  “  X|/2  —  0.019  cm 

Longitudinal  Gamma 

r.-i  -  5.1 

V.-2  “  7 

Beam  Strength  Parameter 

f.  “  0.61 

-  f  r/(4y;2i)  -  5.9  X  10  3 

Transverse  Velocity 

Bru.  i  “  0.135 

Boi  .2  -  0.4  X  10-2 

Critical  Velocity 

Peril.  1  »  0.27 

Peru.  2  “  165  X  10"2 

Output  Wavelength 

X)  “  0.038  cm 

X2  “  2.0  Mm 

Spatial  Growth  Rate 

rmut.i  -  0.37  cm-' 

rmuv.2  “  0.13  cm"1 

Efficiency 

TJ,  -  12% 

y)2  -  0.085% 

Output  Power 

P0i  “  3.6  GW 

P02  ”  25  MW 
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FEL  CONFIGURATION 


OUTPUT  RADIATION 


Figure  1  —  Schematic  of  the  free  electron  laser  model.  The  adiabatic  build-up  of  the  pump  field  is 
not  shown,  and  occurs  to  the  left  of  the  figures  where  the  unmodulated  beam  enters. 


SCHEMATIC  OF  A  TWO  STAGE  FEL 
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Figure  2  —  Schematic  of  the  two-stage  free  electron  laser  concept.  The  electron  beam  enters  at 
left.  Radiation  scattered  at  wavelength  A,  from  the  static  magnetic  pump  of  wavelength  /  in  the 
first  stage  is  reflected  to  act  as  an  electromagnetic  pump  in  the  second  stage.  The  final  scattered  ra¬ 
diation  is  at  wavelength  Aj—  //(*  yf). 
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